Ab Initio Studies of Mechanical, Dynamical, and Thermodynamic Properties of Fe-Pt Alloys

The density functional theory (DFT) framework in the generalized gradient approximation (GGA) was employed to study the mechanical, dynamical, and thermodynamic properties of the ordered bimetallic Fe-Pt alloys with stoichiometric structures Fe3Pt, FePt, and FePt3. These alloys exhibit remarkable magnetic properties, high coercivity, excellent chemical stability, high magnetization, and corrosion resistance, making them potential candidates for application in high-density magnetic storage devices, magnetic recording media, and spintronic devices. The calculations of elastic constants showed that all the considered Fe-Pt alloys satisfy the Born necessary conditions for mechanical stability. Calculations on macroscopic elastic moduli showed that Fe-Pt alloys are ductile and characterized by greater resistance to deformation and volume change under external shearing forces. Furthermore, Fe-Pt alloys exhibit significant anisotropy due to variations in elastic constants and deviation of the universal anisotropy index value from zero. The equiatomic FePt showed dynamical stability, while the others showed softening of soft modes along high symmetry lines in the Brillouin zone. Moreover, from the phonon densities of states, we observed that Fe atomic vibrations are dominant at higher frequencies in Fe-rich compositions, while Pt vibrations are prevalent in Pt-rich.


Introduction
Ferromagnetic Fe-Pt intermetallic alloys are a significant group of solid-state materials, exhibiting intriguing physical and chemical properties that are inherently dependent on their specific chemical compositions.These materials find applications across various industries, including catalysis for chemical reactions, ultrahigh-density magnetic data storage such as hard disk drives and magnetic sensors, spintronic devices, and high-temperature mechanics due to their large uniaxial magnetocrystalline anisotropies (MCA) and high coercivity, amongst other desirable magnetic properties [1][2][3][4].The large MCA presents the potential of decreasing the magnetic particle dimensions to a few nanometers of grain sizes without degrading the thermal stability of the magnetization axis direction.Crystallographically, they stabilize in two Strukturbericht symbols and three space groups: L1 0 (body-centered tetragonal P4/mmm-FePt) and L1 2 (face-centered cubic Pm3m-Fe 3 Pt, Pm3m-FePt 3 and body-centered tetragonal I4/mmm-Fe 3 Pt).The L1 0 phase features an alternating arrangement of Fe and Pt atoms which can undergo a disorder-order phase transition from face-centered cubic (fcc) to tetragonal at temperatures of approximately above 650 • C [5].Müller et al. employed both analytic bond-order potential and latticebased Monte Carlo simulations and reported that the observed disorder-order transition temperature decreases with particle grain size [6,7].Furthermore, Yu et al. reported that the ferromagnetic orientation and dynamical stability break down above the threshold pressure of 96.7 GPa as a result of undergoing spontaneous magnetization [8].The transformation has significant implications for understanding the dynamical properties of Fe-Pt alloys, as it influences the material's magnetic behavior band mechanical response.Furthermore, in L1 0 bimetallic compounds, the difference in atomic size between the two constituent atoms is less than 15%.Thus, the orientation of magnetization and magnitude of the local magnetic moment and stability are significantly influenced by the concentrations and atomic sizes of Fe and Pt [9].These dissimilarities in crystal lattice arrangement and chemical compositions significantly impact the lattice dynamics of these alloys.The face-centered Fe-rich Pm3m-Fe 3 Pt has a similar intermetallic arrangement as the Cu 3 Ag-ordered crystal structure.Temperature-based X-ray diffraction and transmission electron microscopy studies demonstrated that the Pm3m-Fe 3 Pt alloy fully crystallizes at annealing temperatures above 800 °C [10,11].The grain size and lattice parameter increase until this temperature and remain constant thereafter, signifying that crystallization is fully converged and lattice has reached stability.Moreover, the Mössbauer spectrometry and magnetic characterization results showed that the structure conforms to a ferromagnetic state at temperatures below 120 K and paramagnetic state at 300 K and above.The L1 2 cubic Pt-rich Pm3m-FePt 3 is an intermediate alloy that occurs below 300 °C during the formation of the L1 0 FePt [4].Further structural arrangements were discussed in detail in our previous communication [3].Moreover, recently revised phase diagram studies have revealed the existence of other novel ordered FePt 2 and Fe 2 Pt alloys [12].Yu et al. have since performed computational simulations on these alloys to investigate their various thermodynamic, mechanical, and dynamical stabilities [13,14].For the FePt 2 composition, they reported on five different space groups, namely: I4/mmm, P4/nmm, Immm, Cmcm, and P63/mmc, which were predicted to be thermodynamically and mechanically stable, with the I4/mmm being the most thermodynamically stable [13].Furthermore, four space groups of Fe 2 Pt, namely: Cmcm, Immm, I4/mmm, and P3m1 were also reported to be low-energy ferromagnetic stable structures [14].
Phonon dispersion spectra measurements of Fe-Pt alloys to determine their dynamical properties have previously been reported using both experimental and density functional theory (DFT)-based first-principles studies.Pierron-Bohnes et al. measured the vibrational modes of the conventional unit cell of the ordered L1 0 P4/mmm-FePt structure at 300 K using the inelastic neutron scattering method [15].They observed that this alloy is thermodynamically stable since the dispersion relations curves displayed only positive modes of vibrations.Earlier, Noda et al. also investigated and compared the dispersion spectra between the ferromagnetic Pm3m-Fe 3 Pt and paramagnetic Pm3m-FePt 3 cubic L1 2 alloys using the inelastic neutron scattering technique at room temperature [16].They observed that the acoustic modes for both alloys are similar, particularly in the lower-energy region, mainly because their lattice constants differ by only 3%.It has been also determined that the dispersion curves in Fe 3 Pt and FePt 3 are influenced by the force constants associated with their respective nearest neighbors, namely Fe-Fe and Pt-Pt; thus, the nearest neighbor pairs have dominant roles in lattice dynamics.Sternik et al. corroborated the experimental measurements by conducting first-principles computations on the phonon dispersion curves of the 2 × 2 × 2 primitive supercells of P4/mmm-FePt, Pm3m-Fe 3 Pt, and Pm3m-FePt 3 alloys using the direct method as embedded in the PHONON code [17].They reported that the P4/mmm-FePt and Pm3m-FePt 3 alloys show no imaginary modes along high symmetry directions of the Brillouin zone, while there are imaginary frequencies for the transversal acoustic mode on the Pm3m-Fe 3 Pt alloy.Furthermore, their investigation into the impact of magnetic ordering on the dispersion relations of Pm3m-FePt 3 revealed that the frequencies of the antiferromagnetic orientation slightly surpass those of the ferromagnetic orientation.Nevertheless, the dispersion curves exhibit a remarkable similarity.
Investigation of mechanical properties such as elasticity, hardness, shear resistance, ductility, and anisotropy is a crucial prerequisite in solid-state materials before their design and development.DFT calculations of elastic constants have been performed for various binary and ternary alloys and were reported to reveal valuable insights at both ambient and elevated conditions.Phasha et al. investigated the structural and mechanical properties of the ordered Mg-Li binary alloys to predict their phase stability [18].Their findings revealed that there is a correlation between the heats of formation and elastic constants in terms of predicting stability.Phases with negative heats of formation tend to satisfy the Born necessary stability conditions and have a positive tetragonal shear modulus C ′ .Recent DFT studies on various alloys and other compounds have also corroborated these pioneering conclusions.Botha et al. investigated the relationship between the mechanical, dynamical, and thermodynamic properties of Pt-Pd alloys by calculating their elastic constants, phonon dispersion spectra, and heats of formation, respectively [19].Their results indicated that energetically favorable, i.e., with negative heats of formation, compositions satisfy Born stability conditions (mechanical stability) and depict only positive modes of vibrations in the phonon dispersion spectra (dynamical stability).Similar studies were performed on Co-alloyed MnPt alloys, revealing that energetically favorable compositions are mechanically stable [20,21].It was further revealed that the Debye temperature of CsGaSb 2 calculated from the single-crystal elastic constants yields approximately a similar value to the one calculated from the phonon density of states, confirming a robust link between mechanical and dynamical properties.Moreover, the polycrystalline bulk modulus calculated from the elastic constants exhibited exceptional alignment with values derived from alternative equations of state [21].DFT calculations on the elastic constants of Fe-Pt alloys (P4/mmm-FePt, Pm3m-Fe 3 Pt, and Pm3m-FePt 3 ) have also been previously performed [11].It was reported that Pm3m-Fe 3 Pt is mechanically unstable due to negative C ′ and C 44 values, while others satisfied the Born stability conditions.On the contrary, earlier results revealed mechanical stability [7,22].
Despite the significant experimental and computational advances made in these alloys, there remains a notable gap regarding the origins of lattice dynamics and elastic behavior for optimizing their functionality and ensuring their long-term stability.Moreover, little work has been completed on the I4/mmm-Fe 3 Pt alloy, prompting the need to understand and compare the mechanical, dynamical, and thermodynamical properties across the entire spectrum of existing Fe-Pt alloys.Thus, the current communication will put special emphasis on this alloy.These properties are critical in designing and developing new composition-dependent magnetic materials.Additionally, understanding these properties provides valuable insights in predicting the stability of these materials when subjected to industrial conditions.Building from our previous communication [3] where the structural, magnetic, and electronic properties of P4/mmm-FePt, I4/mmm-Fe 3 Pt, Pm3m-Fe 3 Pt, and Pm3m-FePt 3 were reported, the current study presents density functional theory (DFT)-based first-principles computations to establish a corroboration between the mechanical, dynamical, and thermodynamic properties of the same alloys at ambient conditions.Elemental contributions of Fe and Pt atoms on the phonon dispersion spectra are also determined by the partial phonon density of states.Moreover, this work is an extrapolation of the existing theoretical and inelastic neutron scattering measurements of magnetic, mechanical and lattice dynamical properties of Fe-Pt alloys which have reported TA-mode softening on Pm3m-Fe 3 Pt [11,17,23,24].We found that the equiatomic P4/mmm-FePt alloy exhibits dynamical stability, as evidenced by only positive vibrations in the phonon dispersion spectrum.Moreover, the less-studied I4/mmm-Fe 3 Pt alloy shows comparable mechanical stability, ductility, anisotropy, and Debye temperature with the other alloys.Our findings contribute to the growing body of knowledge in the field of alloy research and provide a solid foundation for future investigations.By unravelling the complexities of these ordered alloys, we can unlock their full potential and harness their properties for technological advancements.

Computational Procedure
All calculations were performed using the Cambridge Serial Total Energy Package (CASTEP) [25] simulation code which is incorporated in with the Materials Studio 2020 version suite.CASTEP is a widely used code in the field of computational chemistry and materials science due to its ability to accurately predict the properties and behaviors of materials at the atomic level.Its key strength lies in its ability to perform first-principles calculations based on density functional theory, which accurately describes the electronic structure of materials without the need for empirical parameters, thus allowing the study of properties such as elastic constants, phonon spectra, and thermodynamic properties.However, some calculations are time-consuming and computationally expensive, limiting the size of systems that can be studied and the timescales that can be simulated.We carried out collinear spin polarization density functional theory (DFT) computations within the generalized gradient approximation (GGA) utilizing the Perdew, Burke, and Ernzerhof (PBE) functional [26].All structures were initially relaxed by performing a full geometry optimization while maintaining a fixed basis quality and allowing cell volume to change.The Coulomb interactions between the valence electrons of iron (Fe: 3d 6 4s 2 ) and platinum (Pt: 4f 14 5d 9 6s 1 ) atoms, and their respective pseudo-ionic cores, were accurately modelled using on-the-fly generated (OTFG) ultrasoft pseudopotentials.This approach enables calculations to be conducted with reduced energy cut-offs and ensures that the generated potentials are steady across solid-state and pseudo-atom calculations.Furthermore, it enhances the accurateness and dependability of the results by utilizing the same exchange-correlation functional throughout.A customized energy cut-off of 350 eV was adequate to minimize the total energy of the bulk structure until the difference between two successive low-memory Broyden-Fletcher-Goldfarb-Shanno (BFGS) iterations was within 0.001 eV.The low-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) [27] algorithm was preferred due to its ability to accelerate geometry optimization and accuracy for large systems [28].To sample the wave functions, we employed the Monkhorst-Pack grid parameters of 9 × 9 × 7, 4 × 4 × 4, 6 × 6 × 6, and 8 × 8 × 8 for P4/mmm-FePt, I4/mmm-Fe 3 Pt, Pm3m-FePt 3 , and Pm3m-Fe 3 Pt, respectively.The calculations for all structures were conducted assuming the ferromagnetic (FM) arrangement of the local magnetic moments.Moreover, for the tetragonal P4/mmm-FePt and I4/mmm-Fe 3 Pt, we computed the primitive tetragonal cell as opposed to the conventional cell due to minimal computational resources.As a result, I4/mmm-Fe 3 Pt with lattice parameters a ̸ = c was computed as I4/mmm-Fe 3 Pt with a = c.The phonon dispersion spectra along high symmetry lines and the accompanying phonon density of states were calculated via the finite displacement method which was proven to be highly effective for metallic systems [29,30].The finite displacement technique relies on numerically differentiating forces acting on atoms, which are calculated for multiple unit cells with atomic displacements.Phonon density of states also serves as a requirement for the computation of thermodynamic properties, which permits assignment of temperature factors during analysis.Lastly, the monocrystalline elastic constants were calculated using the stress-strain approach with a maximum strain amplitude of 0.003.The additional calculation criteria for elastic constants involved sustaining total energy of convergence below 2.0 × 10 −6 eV/atom, making sure that the Hellman-Feynman force is maintained under 0.006 eV/Å, and keeping the ionic displacement within 2.0 × 10 −4 .

Mechanical Properties
Before determining the mechanical and dynamical properties of Fe-Pt alloys, structural relaxation was performed to obtain ground-state lattice parameters.To validate the approach utilized, the calculated structural lattice parameters were compared with the existing experimental data.The purpose of this analysis is to ensure the accuracy and reliability of our results.The current GGA-calculated lattice parameters exhibit an impressive agreement of over 98% with the previously reported data, thereby demonstrating the robustness of the approach employed.
Table 1 presents the calculated lattice constants, elastic constants, moduli, Pugh ratio, anisotropy factor, Poisson ratio, Vickers hardness, and Debye temperature for the Fe-Pt alloys.Calculation of elastic constants is crucial in characterizing the mechanical behavior of materials and response to externally induced stress.The Taylor expansion of the total energy of a strained system [31,32] (see Equation ( 1)) was used to calculate the elastic constants (C ij ).
where U(V 0 , 0) is the energy of the unstrained system with equilibrium volume V 0 ; τ i and ξ i are elements in the stress tensor and a factor taking care of the Voigt index, respectively.Cubic and tetragonal crystal systems contain three (c 11 , c 12 , c 44 ) and six (c 11 , c 33 , c 44 , c 66 , c 12 ,c 13 ) independent elastic constants, respectively.Except for C 12 in P4/mmm-FePt, the obtained monocrystalline elastic constants are in agreement with the previous semi-empirical Monte Carlo results obtained using the modified embedded atom method (MEAM) and angular dependent analytic bond-order potential (ABOP) formalism [6,22], further affirming the accuracy of the approach employed.Moreover, our C 12 value is in better agreement with the DFT value of 94 GPa reported by Zotov and Ludwig [11].For cubic crystals (Pm3m-Fe 3 Pt and Pm3m-FePt 3 ) to be deemed mechanically stable, the mandatory Born stability criterion in Equation ( 2) must be fulfilled [33,34].The mechanical stability conditions for tetragonal systems (P4/mmm-FePt and I4/mmm-Fe 3 Pt) are delineated based on Equation (3) [35].We note that the elastic constants for both the cubic and tetragonal systems are satisfied, demonstrating mechanical stability, which corroborates the thermodynamic stability reported previously [3].Additionally, we have noted a positive C ′ value, which further confirms the mechanical stability of the systems.
The elastic constants of the tetragonal systems can further be analyzed by examining their similarities and discrepancies.Significant differences are observed between the calculated C ij values, suggesting that the tetragonal alloys are highly anisotropic.Particularly, C 11 and C 33 are significantly higher than C 44 and C 66 , indicating greater resistance to unidirectional compression compared to shear.Furthermore, there exists a substantial difference between the following pairs of C ij s, C 11 and C 33 ; C 44 and C 66 ; C 12 and C 13 , which has an impact on the behavior of the material.C 11 represents the stiffness of the material in the direction perpendicular to the crystallographic planes, while C 33 represents the stiffness parallel to the crystallographic planes, namely [100]; [010]; and [001].C 11 values are significantly higher than those of C 33 , indicating greater resistance to compression in the a-axis in comparison to the c-axis or greater linear compressibility along the c-axis compared to the a-axis when the material is subjected to external forces.Thus, the interatomic chemical bonds within the (001) plane are more robust compared to the bonding along the crystal direction [001].The relationship C 44 > C 66 , indicates that shear deformation is more prominent when a net stress is imposed along the [100] crystallographic direction of the (010) crystallographic plane, as opposed to when the same stresses are exerted in the [100] direction within the plane (001).This indicates that the alloys exhibit anisotropic behavior with greater resistance to shear deformation in one direction compared to another.Moreover, the significant difference between C 12 and C 13 demonstrates that when stress is imposed along the a-axis, the subsequent strain will be more pronounced along the a-axis than the c-axis.This further highlights the anisotropic response of these alloys under stress, as evidenced by the high A U values.When assessing the level of anisotropy, we have chosen to utilize the universal anisotropy index A U in Equation (4) as postulated by Ranganathan and Ostoja-Starzewski [36].This index is particularly suitable as it considers both shear and bulk contributions, acknowledging the inherent anisotropy present in all single crystals.The magnitude of the deviation of A U from zero is the measure of the degree of anisotropy exhibited by a single crystal, with a value of zero indicating a locally isotropic crystal.
where G V , G R , B V , and B R are the shear and bulk moduli Voigt and Reuss estimates, respectively.We observe that A U decreases with Fe content, thus Pm3m-Fe 3 Pt and I4/mmm-Fe 3 Pt are highly anisotropic over P4/mmm-FePt and Pm3m-FeP 3 , respectively.The calculation of independent elastic constants allowed us to estimate the macroscopic elastic bulk (B), shear (G), and Young's (E) moduli utilizing the Voigt [37] and Reuss [38] estimates.These estimates establish simple and linear relationships between the isotropic bulk and shear moduli of the polycrystalline material.Additionally, Hill demonstrated that the Voigt and Reuss expressions serve as upper and lower bounds, respectively, and proposed an arithmetic average modulus value [39].The bulk modulus signifies the material's ability to withstand compression under applied hydrostatic pressure, while the shear modulus indicates its resistance to deformation from external forces at a constant volume.Young's modulus, on the other hand, measures the material's stiffness by determining the ratio of vertical linear stress to linear strain.
The current DFT calculations revealed that Fe-Pt alloys have intrinsic mechanical hardness, stiffness, and shear resistance, i.e., are highly resistant to changes from external mechanical stress due to relatively large magnitudes of B, G, and E. The shear modulus is less than the bulk in all the alloys, suggesting they exhibit greater resistance to volumetric change than shear deformation, and that the shear modulus is the parameter limiting stability.Pm3m-FePt 3 possesses the highest bulk modulus value, suggesting the greatest mechanical hardness and compression resistance over P4/mmm-FePt, Pm3m-Fe 3 Pt, and I4/mmm-Fe 3 Pt, respectively.This is consistent with a previous DFT study on Pt/Pd alloys, which showed that the bulk modulus is higher in Pt-rich compositions [19].The Pm3m-Fe 3 Pt system possesses the lowest shear modulus value of 46.379 GPa, indicating the greatest susceptibility to shear deformation compared to I4/mmm-Fe 3 Pt, Pm3m-FePt 3 , and P4/mmm-FePt, respectively.The Young's modulus suggests that P4/mmm-FePt adopts highest stiffness over, Pm3m-FePt 3 , Pm3m-Fe 3 Pt, and I4/mmm-Fe 3 Pt, respectively.
To assess ductility and brittles, we computed the Pugh shear to the bulk modulus (K) [40] and Poisson (υ) [41] ration of solid-state materials.Pugh postulated that a material is considered ductile if K is less than the critical value 0.5 and brittle when greater than 0.5.We noticed that K values are less than 0.5 for all Fe-Pt alloys, demonstrating ductility.Ductile materials exhibit greater resistance to thermal shock and are easier to machine.Moreover, we noted that the Fe-rich Pm3m-Fe 3 Pt and I4/mmm-Fe 3 Pt systems possess lower K values, due to their low shear modulus.The Poisson's ratio associated with volume changes in directions perpendicular to the applied tension during deformation typically ranges from −1 to 0.5.An indication of ductility is present when the value exceeds 0.26, otherwise the material is considered brittle.The Poisson ratio for these alloys is greater than 0.26, further confirming ductility.Interestingly, the less-studied I4/mmm-Fe 3 Pt shows excellent ductility over Pm3m-FePt 3 and P4/mmm-FePt, respectively.Furthermore, the Poisson ratio was used to analyze bonding type in Fe-Pt alloys.Crystals are considered covalent if υ is approximately 0.1, ionic if υ = 0.25, and metallic if υ > 0.33 [42,43].We note that υ values are greater than 0.33 for Pm3m-Fe 3 Pt, I4/mmm-Fe 3 Pt, and Pm3m-FePt 3 , implying metallic bonding.On the other hand, P4/mmm-FePt shows a ratio less than 0.33 but greater than 0.25, indicating half-metallic behavior.This is in good agreement with the density of states predictions reported in our previous communication [3].To gain insight into the thermal properties of Fe-Pt alloys, we determined the Debye temperature (θ D ).This property represents the point at which the atoms within a solid material cease to vibrate independently, but instead begin to move in a synchronized manner.It signifies the highest normal mode of vibration and serves to establish a correlation between elastic constants and phonon dispersions, specific heat, thermal expansion, and thermal conductivity.Equation ( 8) was utilized in the computation of the Debye temperature [44].
where n, ν, and k B are the atom concentration, phase velocity, and Boltzmann constant, respectively.Our DFT calculations show relatively huge magnitudes of Debye temperature (~300 K), indicating high vibrational modes in the phonon dispersion spectra (see Figure 1) and consequently high thermal conductivity.Laureti et al. employed the Correlated Debye Model (CDMT) to measure the Debye temperature of P4/mmm-FePt and reported θ D = 340 K [4], which is consistent with our results (330.79K). Materials with higher Debye temperatures also have stiffer and more rigid lattices, while those with lower Debye temperatures exhibit softer and more flexible structures [45,46].To estimate the Vickers hardness of Fe-Pt alloys, we employed the microscopic bond resistance model proposed by Tian et al. [47] as shown in Equation ( 9).
where k is the G/B ratio and G is the shear modulus.Our findings show that the metallic Fe-Pt alloys relatively show low hardness.In principle, this is expected for ductile materials with low (<0.5)G/B values, as proposed by Chen et al. [48].
lower Debye temperatures exhibit softer and more flexible structures [45,46].To estimate the Vickers hardness of Fe-Pt alloys, we employed the microscopic bond resistance model proposed by Tian et al. [47] as shown in Equation ( 9). = 0.92 . ., (9) where  is the   ⁄ ratio and  is the shear modulus.Our findings show that the metallic Fe-Pt alloys relatively show low hardness.In principle, this is expected for ductile materials with low (<0.5)  ⁄ values, as proposed by Chen et al. [48].

Phonon Dispersion Curves
The phonon dispersion relations along high symmetry lines in the Brillouin zone were calculated to determine the lattice dynamic stability of Fe-Pt alloys and are presented in Figure 1.Phonon dispersion frequencies are calculated as the analysis of forces associated with a systematic set of atomic displacements from equilibrium positions in a crystal system.Materials are considered stable if there are no negative frequencies (soft modes) along high symmetry lines.We observe that the ordered L10 4/-FePt displays only positive phonon vibrational modes along the high symmetry directions, indicating that all the eigenvalues are real and the alloy is dynamically stable.This finding aligns well with previous theoretical and experimental research [15,17,49].The other compositions, in particular 3 -FePt3, show imaginary frequencies of vibration along high symmetry lines in the spectrum, reflecting deviation dynamical instability and possible structural deformation which results in lowering of the symmetry [50].The anomalous behavior shows that the frequencies of vibration decrease with the increasing wave vector as opposed to increasing.This may be attributed to the presence of complex interactions between the Fe and Pt atoms in the lattice, leading to the formation of localized vibrations that propagate through the material in a non-traditional manner.Dynamical instability is

Phonon Dispersion Curves
The phonon dispersion relations along high symmetry lines in the Brillouin zone were calculated to determine the lattice dynamic stability of Fe-Pt alloys and are presented in Figure 1.Phonon dispersion frequencies are calculated as the analysis of forces associated with a systematic set of atomic displacements from equilibrium positions in a crystal system.Materials are considered stable if there are no negative frequencies (soft modes) along high symmetry lines.We observe that the ordered L1 0 P4/mmm-FePt displays only positive phonon vibrational modes along the high symmetry directions, indicating that all the eigenvalues are real and the alloy is dynamically stable.This finding aligns well with previous theoretical and experimental research [15,17,49].The other compositions, in particular Pm3m-FePt 3 , show imaginary frequencies of vibration along high symmetry lines in the spectrum, reflecting deviation dynamical instability and possible structural deformation which results in lowering of the symmetry [50].The anomalous behavior shows that the frequencies of vibration decrease with the increasing wave vector as opposed to increasing.This may be attributed to the presence of complex interactions between the Fe and Pt atoms in the lattice, leading to the formation of localized vibrations that propagate through the material in a non-traditional manner.Dynamical instability is not pronounced for Pm3m-Fe 3 Pt and I4/mmm-Fe 3 Pt since the negative vibrations are not along the center of the Brillouin zone G.
The modes of vibration for Pm3m-Fe 3 Pt and I4/mmm-Fe 3 Pt as depicted in Figure 1a,b share a related frequency range and similarities along G, which can be attributed to their identical chemical composition and stoichiometry.Moreover, we note that the separation between the acoustic and optical modes is not enough to create a gap in the dispersion relations of all the alloys, suggesting seamless continuous energy transfer between the phonon modes.The disparities in mass among the alloys are evident through the observation that the frequencies of vibrations (optical modes) are higher in Fe-rich (with smaller mass) alloys and lower in Pt-rich (with larger mass) alloys.

Phonon Density of States
To analyze the distinctive elemental vibrational contributions, we conducted calculations of the phonon density of states (DOS) as illustrated in Figure 2. The spectra can be categorized into two frequency regions: one where Fe vibrations predominate and another where Pt vibrations are more prominent.In Fe-rich compositions, Fe vibrations are prevalent at higher frequencies, whereas Pt-rich compositions exhibit dominance of Pt vibrations.Additionally, a significant level of anisotropy in lattice dynamics is apparent in all structures, particularly in P4/mmm-FePt (Figure 2c).The lowest frequency band (1.9-3.6 THz) comprises vibrations of the heavier platinum atoms, while frequencies above 4 THz are attributed to Fe vibrations.These findings align well with similar DFT calculations reported in previous studies [17,50].In the Fe-rich Pm3m-Fe 3 Pt and I4/mm-Fe 3 Pt compositions, Pt movement prevails up to 4 THz with minimal contribution from Fe, indicating that heavier Pt movement is responsible for the presence of imaginary soft modes leading to dynamical instability.Conversely, in the Pm3m-FePt 3 composition, Fe vibrations dominate in the −5 to 2 THz range, while Pt is more prominent in the 4 to 6 THz region.Therefore, the negative frequencies observed in the phonon dispersion curves of this alloy stem from the movement of multiple Fe atoms.Moreover, it should be noted that the elemental partial DOSs of Fe and Pt were sampled for only one atom each.Hence, there are significant differences between the sum (red) and atomic plots (blue and green) in the non-equiatomic compositions since other atoms are not included.

Thermodynamic Properties
Thermodynamic (temperature-dependent) quantities have been evaluated from the phonon density of states (DOS) in the temperature range of 0 to 1000 K based on the relations developed by Baroni et al. as listed in Table 2 [51].All quantities were calculated at ground state, that is, geometry optimization was fully converged and all the phonon eigenfrequencies were real and non-negative.

Thermodynamic Properties
Thermodynamic (temperature-dependent) quantities have been evaluated from the phonon density of states (DOS) in the temperature range of 0 to 1000 K based on the relations developed by Baroni et al. as listed in Table 2 [51].All quantities were calculated at ground state, that is, geometry optimization was fully converged and all the phonon eigenfrequencies were real and non-negative.Enthalpy, free energy, and entropy are fundamental thermodynamic properties that govern the behavior of solid-state materials.In the context of solid-state materials, enthalpy plays a crucial role in understanding phase transitions, such as melting and crystallization.For example, the enthalpy change associated with the melting of a solid reflects the energy required to overcome intermolecular forces and disrupt the crystalline structure, leading to a transition from a solid to a liquid state.Free energy is a key parameter in predicting the stability and equilibrium of different crystal structures.By comparing the free energy of different crystal phases, we can determine the phase that is thermodynamically favored under specific conditions.Entropy is closely related to the organization of atoms and molecules in a crystal lattice.Higher entropy is associated with increased disorder and greater freedom of movement for particles, whereas lower entropy corresponds to a more ordered and structured arrangement which is associated with favorable thermal conductivity and resistance to thermal expansion and deformation.
These three thermodynamic potentials as a function of temperature are plotted in Figure 3. From our graphical representations, we observed that entropy and enthalpy rise with temperature, while the free energy decreases.The I4/mmm-Fe 3 Pt alloy shows the largest entropy values, alluding to the high mobility of atoms within the lattice structure, while the P4/mmm-FePt has the lowest, suggesting a more compact structural arrangement and thermal stability and conductivity.This is consistent with the phonon dispersion spectra predictions where P4/mmm-FePt showed dynamical stability.The free energies are negative and continue to decrease with temperature, implying thermodynamic stability [49], which corroborates the heats of formation in the previous communication [3]. Figure 4a shows the variation in isobaric heat capacity ( C v ) as a function of temperature (0-1000 K).Heat capacities are critical in determining the thermal behavior of solid-state materials and are strongly dependent on intermolecular bonds.Materials with strong intermolecular bonds tend to have higher heat capacities as more energy is required to disrupt these bonds and increase the temperature.As the temperature increases, the heat capacities of Fe-Pt alloys show uneven characteristics.At temperatures below 200 K, C v increases exponentially and gradually scales off toward the Dulong-Petit boundary at higher temperatures.This behavior is common in most solid-state materials [19,21,52].In the temperature range, I4/mmm-Fe 3 Pt (~46 cal/cell•K) shows the highest C v value, followed by Fm3m-Fe 3 Pt (~23 cal/cell•K), Fm3m-FePt 3 (~17 cal/cell•K), and P4/mmm-FePt (~12 cal/cell•K), respectively.Interestingly, in temperatures below 80 K, the heat capacities are consistent with the bond length (see Table 1).An alloy with the shortest bond length (I4/mmm-Fe 3 Pt) exhibits the highest C v value, while an alloy with the longest bond length (P4/mmm-FePt) has the lowest value, further confirming a compact lattice structure on the latter.

Conclusions
The current study has successfully conducted ab initio simulations on the bimetal Fe-Pt alloys to gain insights into the intrinsic factors underlying their mechanical, dyna ical, and thermodynamic behavior, as well as to derive their stability trends.Moreov we have established a link between the lattice dynamic and thermodynamic properti All Fe-Pt alloys were predicted to be mechanically stable since they satisfied the Born s bility conditions for cubic and tetragonal crystal lattices.Moreover, Fe-Pt alloys are ch acterized by the presence of sizeable anisotropy rising from the discrepancies in elas constants and deviation of the universal anisotropy index from zero, which may give r to a type of anisotropy called magnetocrystalline anisotropy.The prediction of mechani stability and anisotropy is in good agreement with the thermodynamic stability a

Conclusions
The current study has successfully conducted ab initio simulations on the bimeta Fe-Pt alloys to gain insights into the intrinsic factors underlying their mechanical, dyna ical, and thermodynamic behavior, as well as to derive their stability trends.Moreov we have established a link between the lattice dynamic and thermodynamic propert All Fe-Pt alloys were predicted to be mechanically stable since they satisfied the Born s bility conditions for cubic and tetragonal crystal lattices.Moreover, Fe-Pt alloys are ch acterized by the presence of sizeable anisotropy rising from the discrepancies in ela constants and deviation of the universal anisotropy index from zero, which may give r to a type of anisotropy called magnetocrystalline anisotropy.The prediction of mechan stability and anisotropy is in good agreement with the thermodynamic stability a The variation in the Debye temperature ( θ D ) as a function of temperature is shown in Figure 4b.Our results show a very strong linear increment relation between θ D and T for Pm3m-FePt 3 and a moderate linear increase for I4/mmm-Fe 3 Pt and Pm3m-Fe 3 Pt.Interestingly, P4/mmm-FePt shows a nearly constant θ D vs. T relationship with a small standard deviation within 14, suggesting that the modes of vibration are compact and thermal stability is achieved for P4/mmm-FePt, which is in good agreement with the phonon dispersion spectra and entropy predictions.The Debye temperature decreases with Pt content, which is in good agreement with related DFT studies [19,53].Moreover, θ D of Pm3m-FePt 3 is highly affected by temperature, indicating the increase in atomic mobility and vibrations, which is consistent with the negative phonon modes of vibration.

Conclusions
The current study has successfully conducted ab initio simulations on the bimetallic Fe-Pt alloys to gain insights into the intrinsic factors underlying their mechanical, dynamical, and thermodynamic behavior, as well as to derive their stability trends.Moreover, we have established a link between the lattice dynamic and thermodynamic properties.All Fe-Pt alloys were predicted to be mechanically stable since they satisfied the Born stability conditions for cubic and tetragonal crystal lattices.Moreover, Fe-Pt alloys are characterized by the presence of sizeable anisotropy rising from the discrepancies in elastic constants and deviation of the universal anisotropy index from zero, which may give rise to a type of anisotropy called magnetocrystalline anisotropy.The prediction of mechanical stability and anisotropy is in good agreement with the thermodynamic stability and sizeable magnetocrystalline anisotropy energies reported in our previous communication.Computation of phonon dispersion spectra showed that the equiatomic L1 0 P4/mmm-FePt is dynamically stable at ambient conditions, while the other compositions show soft modes along high-symmetry directions of the Brillouin zone.This study revealed that the negative vibrational modes predominantly emanate from Fe atoms in the highly dynamically unstable Pm3m-FePt 3 alloy, while dynamical instability is not pronounced in I4/mmm-Fe 3 Pt since the negative vibrations are not along the center of the Brillouin zone.Moreover, I4/mmm-Fe 3 Pt shows excellent mechanical stability, ductility, and thermal properties.Furthermore, the temperature versus entropy plots indicated that dynamically unstable alloys exhibit high atomic mobility as entropy exponentially increases with temperature, leading to the disordering of the lattice structure.This was corroborated by the heat capacity and Debye temperature plots.The free energies were predicted to be negative and to continue to decrease with temperature, implying thermodynamic stability in all the alloys.In addition to the current conclusions, it is of interest to perform ternary alloying of FePt with other transition metals such as Mn, Co, and Ni to determine any potential enhancement of the desired properties.

Table 1 .
Calculated lattice constants, elastic constants, and moduli for the four ordered bimetallic Fe-Pt alloys.Available experimental data are provided in parentheses.